Pinning complex Networks by A Single Controller 
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Abstract In this short paper, we point out that a single 
local stability controller can pin a linear or nonlinear coupled 
complex network to a specified solution (or an equilibrium) 
of the coupled complex network. A rigorous mathematical 
proof is given, too. 

Keywords: Complex networks, chaos synchronization, ex- 
ponential stability 

I. Introduction 

Many natural and man-made systems, such as neural 
systems, social systems, WWW, food webs, electrical power 
grids, etc., can all be described by graphs. In such a graph, 
every node represents an individual element of the system, 
while edges represent relations between nodes. These graphs 
are called complex networks. For decades, complex networks 
have been focused on by scientists from various fields, for 
instance, sociology, biology, mathematics and physics,etc. 

Linearly Coupled Ordinary Differential Equations 
(LCODEs) are a large class of dynamical systems with 
continuous time and state, as well as discrete space to 
describe coupled oscillators. 

In general, the LCODEs can be described as follows: 



X 



dx (t) m 



»(*)]. 



1, • • • , m 



where Xi(t) — [xj(t), ■ ■ ■ , a;"(i)] T e R n is the state variable 
of the i — th node, t £ [0, +oo) is a continuous time, / : 
R x [0, +oo) — > R n is continuous map, ay > for i.j = 
L ■ • ■ , m, i ^ j. 
■ Letting an = 



rewritten as follows 

dxi(t) 



o-ij, the equation above can be 

j=ij¥* 



dt 



f(xi(t),t) + }^a ij Xj(t) i = 1, • ■ • , to 

3=1 



(1) 



where A = (ctij) £ n my - m [ s the coupling matrix with zero- 
sum rows and <2y > 0, for i ^ j, which is determined by the 
topological structure of the LCODEs. 

Today, in science and technology , "synchronization" ap- 
pears in a wide range of real systems. In theoretical fields, 
there are various kinds of concepts of synchronization. For 
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example, phase synchronization, imperfect synchronization, 
lag synchronization, and almost synchronization, etc. In this 
paper, the concept of complete synchronization is considered, 
which can be mathematically defined as follows: If a coupled 
system composed of m sub-systems 



fifat, ■ ■ 



,x m ,t), for i = l,---,m 
■Xi(t)\\ = 0, for alii, j = 1,- 



, to, then 



satisfies lim 

t— >oc 

the coupled system is said to be completely synchronized, for 
simplicity, synchronized. 

The synchronization of the LCODEs has attracted much 
attention of researchers in different fields. In [1], [2], the 
local stability of the synchronization manifold was studied 
via linearization method named "master stability function". In 
[3], a distance was defined between the state of the coupled 
system and synchronization manifold. Using this distance, a 
methodology was proposed to discuss global convergence for 
complete regular coupling configuration (also see [4]). 

In the paper [5], based on the geometric structure of the syn- 
chronization manifold, a thorough theoretical analysis of the 
synchronization of LCODEs is given. Moreover, it was pointed 
out that the synchronization of the coupled system only means 
that lim ||x 7 -(t) — £j(t)|| = 0, for all i,j — l,---,m. It 

t— >oo 

does not mean that lim Xj(t) — s(t) for some s(t) satisfying 

t— >oo 

s(t) = f(s(t)). Even if all the differences ||xj(0) - s(0)|| are 
very small, it still hold that lim \\xj(t) - s(t)\\ ^ 0. 

t— *oo 

However, in some cases, it is needed that all the states 
of the the nodes synchronize to some special solution s(t) 
of the system s(i) = f(s(t)) or some equilibrium point s* 
satisfying /(a*) = 0. Moreover, it often happens that due to 
some reasons, some nodes can not receive controller, thus, it 
is natural to raise the question: can we realize synchronize all 
states Xi(t) to s(t) only add control to some nodes. 

In the paper [6], the authors discussed how to pin a complex 
dynamical networks, where the coupling matrix is symmetric, 
to its equilibrium by adding some controllers. They pointed 
out that significantly less (but do not specify how many) local 
stable controllers can pin a complex dynamical networks to 
its equilibrium. Here, we rigorously prove that only a single 
controller can achieve this goal. 

II. Pin a complex network with linear coupling 
Suppose that the complex network is 

dx l {t) 



dt 



f (xj (i),t)+cyi ajj Xj (t) i = l,---,m (2) 

3=1 



where Xi G R n , a,tj > 0, i / j and 



0, for i 



3=1 



1, 2, • • • , to. s(t) is a solution of the uncoupled system 

m = f(s(t)) o) 



We prove that if e > and c is large enough. The following 
coupled network with a single controller 



dxi(t) 
dt 

dxjjt) 
dt 



= f(xi{t),t) +c^2aijXj(t) - ce{xi{t) - s(t)), 



f{xi{t),t) + c^2 i a i jXj{t), i = 2,---,m (4) 



can pin the complex dynamical network (1) to s(t). 

Denote Sxi(t) — xt(t) — s(t), then the system (1) can be 
rewritten as: 

= f{x t {t),t) /( S (t)) +cf> J fa J (t), 



i = !,■■■ ,m (5) 



and the network with a single controller @ is written as 

dSxi(t) 



dt 



f(xi(t),t) - f{s(t)) +cJ2^Sx,(t), 

i = 1, ■ • • , m (6) 



where 5n = an — e, e > and a,y = otherwise. 
At first, we prove the following simple proposition. 
Proposition 1 If A = (cty)Fj=i * s a matr i x satisfying 

ra 

a > 0, if i j, and a,j = 0, for i = 1, 2, • ■ • , m. 



an = a 



i=i 



Then, all eigenvalues of the matrix 

/ an - e a 12 

0-21 0,22 



A 



\ a m i a m 2 



aim \ 

0-2m 
Oram J 



are negative. 

Proof Suppose that A is an eigenvalue of A, v = 
[vi, ■ ■ ■ , v m ] T is the corresponding eigenvector, and \vk\ = 
maXj—x t ... tm \vj\. It is clear that if v is an eigenvector, then 
— v is also an eigenvector, thus, without loss of generality, we 
can assume that Vk > and Vk = maxj = i i ... im \vj\. 

If k = 1. Then 



ZVi 



< -evi + dylvjl < — evi < 

which means A < 0. 
Instead, if k > 1. Then, 

m m 

AWfe = },akjVj < dkkVk + ^akjlvjl < 

which means A < 0. If A = 0, then, w = [w/c • • • , ffc] T . 
However, this is impossible. For 

n i 

^ 5ij«j = -ev k < 
Therefore, A < 0. The proposition is proved. 



A. Pin a complex network with linear coupling and symmetric 
coupling matrix 

In this section, we investigate the complex networks, where 
the coupling is linear and the coupling matrix is symmetric. 

With the simple proposition given above, we prove two the- 
orems. Theorem 1 addresses local synchronization. Theorem 
2 addresses global synchronization. 

Let Sx(t) = [Sxx(t), ■ ■ ■ , Sx m (t)]. Differentiating 
along s(t) gives 



dSx{t) 
dt 



= D(f(s(t)))Sx(t) + cox(t)! 1 



(7) 



Let A T = WJW T be the eigenvalue decomposition of A, 
where J — diag{X%,- ■ ■ , X m }, > Ai > ■•• > A m , and 
Sy(t) = Sx(t)W. Then we have 



dSy k (t) 
dt 



[Df(s(t))+cX k I]Sy k (t) 



(8) 



Theorem 1. Let Hi(t), i — 1,- ■ ■ ,m, are the eigen- 
values of the matrix \(Df(s{t)) + Df T (s(t))), /j,(t) = 
maxi = i.... !m fii(t). If n(t) < — cAi for all t > 0. Then, 
the coupled system with a controller (4) can be locally 
exponentially synchronized to s(t). 

Proof: It is easy to see that 



l d{SyJ(t)6y k (t)} 
2 dt 



= SyJ(t)[Df(s(t)) + cX k I]5y k (t) 



= 6yJ(t)[-(Df( S (t)) + Df T (s(t))) + cX k I]5y k (t) (9) 



Under < — cAi for all t > 0, we have 



l d{SyJ(t)Sy k (t)} 
2 dt 



<SyJ(jtMt)-cXi]Sy k (t) (10) 



which means that 6y] (t)Sy k (t) = 0(e^W- cAl )*). Theorem 
1 is proved. 

Theorem 2. Suppose > Ai > • ■ • > X m are the 
eigenvalues of A. If there are positive diagonal matrices P = 
diag{pi, ■ ■ ■ ,p n }, A = diag{A%, ■ ■ ■ , A„} and a constant 
i] > 0, such that 

(x~y) T P{f(x 1 t)-Ax-f{y 7 t) + Ay) < - v (x-y) T (x-y) 

(11) 

and Afc + cAi < for k = 1, ■ • ■ ,n. Then, the controlled 
system (4) is globally exponentially synchronized to s(t). 

Proof: Define a Lyapunov function as 



V(t) = -Y / Sx l (t) T PSx l (t) 



Denote Sx k (t) = [5x\{t), • • • , Sx k n (t)] T . Then, we have 

d ^ = jTSx i (t) T Pf(x i (t),t)5x j (t) 

i—1 

m 

- E fe i(*) Tp /( s (*)) fa 3-(*) 

i=l 

m rn 

+ Sxj(t) T Pc dijSxj(t) 

i=l 3=1 



i=l 



P(/(ii(i),i) -/(«(*))) 

+ ^fa,(() T F c^a lj fc j (()4-Afe 1 (i) 
»=i L j=1 

m 
i=l 

m r m 

»=i L j- =1 

m 
i=l 

n 



matrices P = diag{pi, ■ ■ ■ ,p n }, A = diag{Ai, • • • , A„} and 
a constant 77 > 0, such that 

(x- 2 /) T P(/(x,i)-A 2 ;-/(y,i) + A2;) < -r;( a; -y) T (x- 2 /) 

(13) 

and A/. + acAi < for k = 1, ■ • • ,n. Then, the controlled 
system (II 21 is globally exponentially synchronized to s(t). 

Proof Along with g(x t ) = [g(xj), ■ ■ ■ , g(x™)} T , i = 
l,---,m, we denote g{x k ) = \g(x\), ■ ■ ■ , g(x k n )] T , k = 
1, • • • , n. 

We use the same Lyapunov function 



V(t) = -Y / Sx l (t) T PSx l (t) 



In this case, we have 

dV ^ -±SxJ { t)P dSx ^ 
»=i 



/(**(<),*) - /(a(t),i) 



cE^S^C*)) 
3=1 J 

- C£ [xi(t) - s(t)] T P[ ff ( 2; i(t)) - g(s(t))} 



k=i 

Because cA + A&J < 0, we have 

dV(t) 

^ —'I 

»=i 

Therefore 



< 



< -7 ? Vfe l (i) T (5a; l (t) < . 
at f—f mini Pi 



= 0(e~^T¥-) 

Theorem 2 is proved completely. 

Remark 1 It is clear that if c is large enough, then the 
coupled network with a single controller can pin the complex 
network to a solution s(t) of the uncoupled system. 

B. Pin a complex network with nonlinear coupling 

In this section, we discuss how to pin a complex network 
with nonlinear coupling. In this case, the coupled system (0J 
takes the following form 

=f(xi(t) ) t)+cY / a lj9 (x j (t)) 



+c^a l] g(x J {t)) 
3=1 J 
-ce[xi(i) - s(t)] T Pa[xi{t) - s(t)] 

By the property of the matrix A, it is easy to verify that for 

u = [ui, • • • , u m ] T , v = [«!,••-, w,„] T , 

m 

U T Au = E u i a tj v J = E fl y ( W « _ _ V j) (I 4 ) 

i,j=l i>i 

combining with 9 ^~ 9 ^ > a > 0, we have 



E<v<-, ;/i/'E"'^-0;'M 



3=1 



= e»^(*) t ^(^(*)) 



3=1 



-cefoMt))-^*))), 



= f(x i (t),t)+c'y]aijg(xj(i)) i = 2, • • • , m (12) 



fc=i 
n 

= Ep^E^-^'w-^w)^^))-^^))) 

/c=l ,;">2 

n 

< -sEwE^-^w - **(*))(*.•(*) - **(*)) 

fc=l 3>i 



3=1 



where g(xi(t)) = [ff(xj(t)) ( • • • , g(x?(t))] T and 3 (.t) is a 
monotone increasing function. 

In the following, we will prove that the coupled complex 
network with a single controller d!2t can pin it to a specified 
solution s(t) = f(s(t),t), too. In particular, we prove 

Theorem 3. Suppose > Ai > • • ■ > A m are the eigen- 
values of A, 9< " u lZv > a > 0. If there are positive diagonal 



k=l 



Therefore, 



dV(t) 
dt 



< 



Y,PkSx k (t) T [A k I m + caA]6x k (t) 



fe=i 



Because A^ + acXi < 0, we have 



dt 



mm. Pi 



and 



2*_ 



V(t) =0(e mi -^ 
Theorem 3 is proved. 

C. Pin a complex network with asymmetric coupling matrix 

In practice, indirect graphs are small part. Most of the 
graphs are direct graphs. It means the coupling matrix is asym- 
metric. Therefore, we must investigate pining the complex 
networks, in which the coupling matrix is asymmetric. This is 
the issue we investigate in this section. 
Consider the system 



dx\(t) 
dt ~ 

dxj(t) 
dt 



f(xi(t),t) +cJ~]aijXj(t) - ce(xi(t) - s(t)), 



3 = 1 



f(xi(t),t) +c^2aijXj(t), i = 2,---,m (15) 
where the coupling matrix A is asymmetric but satisfies zero 

Em ri 
j=l a ij = °- 

Let [£i, • • • , £, m ] T be the left eigenvalue of the matrix A. It 
is well know that if A is irreducible, all £j > 0, i = 1, • • • , m. 

Define E = diag[£i, • • • , £ m ]. It is easy to verify that EA + 
A T E is a symmetric matrix and zero row sum. Therefore, by 
the Proposition 1, EA + A T E < 0. 

Theorem 4. Suppose that A is irreducible, > p 1 > 
■ ■ ■ > /U TO are the eigenvalues of (EA + A T E)/2. If there 
are positive diagonal matrices P = diag{pi, ■ ■ ■ ,p n }, A = 
diag{A\, • • • , A„} and a constant rj > 0, such that 

(x-y) T P(f(x,t)-Ax-f(y,t) + Ay) < - v (x - y) T (x - y) 

(16) 

and At max £j + cpi < for k — 1, • • • , n. Then, the 

i— l,"',m 

controlled system (4) is globally exponentially synchronized 
to s(i). 

Proof: In this case, define a new Lyapunov function as 

m 
i=l 

Differentiating V(Sx) f(x,t) £ QUAD(A,P), we have 



i=l 



f( Xi (t),t) - f(s(t),t) 



i=l 



= X-' lr '' :/:/ ' 



t=l 



(/(^(t),t)-Ax l (t))-(/( s (t),i) 



As(i)) +cy^aijfaj(t) + Afei(t) 

m 

m r m 

+ ^2tiSxf(t)P cY^a ij 6x j (t) + A5xi(t) 



< 



3=1 



< -7]V(6x) + Y,P]^ T {t)E(cA + A k E m )5x j {t) 

3 = 1 



v(t) l^A . T/ . 

-r] . — + - > p 5x lT {t) 
mm lPl 2^^ 



3 = 1 



E(cA + A k E m ) 



x5S: 3 (t) < -r?-^ 



mini Pi 

For under the assumption Aj + cp\ < 0, we have 



E(cA + A k E m ) 
Therefore 



is negative definite for k — 1, • • • , n. 



V(i)=0(e 

Theorem 4 is proved completely. 

In the following, we remove assumption that A is irre- 
ducible. In this case, we assume 



A 



An 
A21 A 2 2 

A n i A, 



pi jL ip2 ' ' ' 



(17) 



where A qq G R m <i' m i t q — 1, • • • are irreducible. And, for 
each q, there exists q > k such that A qk ^ 0. 

It is easy to see that if we add a single controller —e(xi (t) — 
s(t) to the node Xi(t). Then, by previous arguments, we 
conclude that the subsystem 



dx\(t) 
dt 



f(x 1 (t),t) + cy^cujXjjt) - ce(xi(t) - s(t)), 



3=1 



dxj(t) 
dt 



f(x i (t),i) + c'y]aijXj(t), i = 2,---,m (18) 
3=1 



pins xi(t), ■ ■ ■ ,x mi to s(t). 



dxjjt) 
dt 



= f(xi(t),t) + c'y]a i jXj(t), i = mi + 1, • • • ,m 2 
3=1 



can be written as 

d5xi(t) 



dt 



= f( Xi {t),t)- f(s(t),t)+c aij^jit) + 0(e 

j=m 1 +l 



Because A 2 i ^ 0. Then, in A22, there exists at least one 
row i 2 , such that 



a 42«2 > y ' a i2j X j(t) 

j=m 1 + l 



(19) 
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Fig. 1. Chaotic Behavior of the Chua's circuit. 

Fig. 2. The coupled system 1211 is synchronized without controller. 



Therefore, similar to the proposition, we conclude that all 
eigenvalues of A22 are negative. By the similar arguments in 
the proof of theorem 4, we can pin Xi(t), i = mi + 1, • • • , rri2, 
to s(t). 

By induction, we prove that if we add a controller to the 
master subsystem corresponding to the sub-matrix An, then 
we can ping the complex network to s(t) even if the coupling 
matrix in reducible. 

III. Simulation 

In this section, we give some numerical examples to 
verify the theorem given in previous section. 
We consider the Chua's circuit 



= k[x 2 - h{xx)\ 

dt 

»■'•:< 

dt 



xi - x 2 + X3 

-lx 2 



(20) 



2™ _ _3 
7 x 



au x+ 1| - \x- 1|], k = 9 and I = Hf. 



where h(x) - _.. 
With these parameters, the system has a double-scroll chaotic 
attractor, as shown in Figure 1. 

As the coupled system, we consider three globally coupled 
Chua's circuits 

3 

f( Xi {t),t)+cJ2aijXj(t) (21) 



dx*(t) 
dt 



where x l (t) = (x[(t), x\{t), x\(t)) T S R 3 , i = 1,2,3, and 
/(•) is defined as in d20"l . 



A. The coupling matrix is symmetric 

Suppose that the coupling matrix without the controller is 

-5.1 5.0 0.1 
A = I 5.0 -11.0 6.0 
0.1 6.0 -6.1 

As for the coupling strength, we pick c = 10. 
Direct calculation shows that 

{x-y) T P(f{x,t)-Ax-f(y 7 t) + Ay) < - V (x-y) T {x-y) 

(22) 

where P = I 3 , A = 10J 3 , 77 = 0.6218. 




Fig. 3. The solution of coupled system 1211 does not converge to s(t) without 
controller. 



Let the initial values be xi(0) = (40.1, 20.2, 30.3) T , 
a; 2 (0) = (20.4, 30.5, 10.6) T , 2:3(0) = (60.7, 40.8, 50.9) T . 
s(t) = is a solution of the uncoupled system. 

Denote x^(t) — (xi(t) + X2(t) + X3(i))/3. We use the 

V 3 _ \\xi(t)-xe(t)\\ 

quantity =H to measure synchronization capa- 

||a:j(0)-x 4 (0)|| 

bility. Figure 2 indicates that the coupled system i2l\ can reach 
synchronization without controller. We also use the quantity 

V 3 _ ||x<(i)-s(i)|| , . 

i^j- 1 to measure if the specified solution sit) of 

the uncoupled system is stable. Figure 3 indicates that s(t) is 
unstable for the coupled system (12 1 i without controller. 

Now, we add a single controller to the first node of the 
coupled system with e = 4.9 



dxj(t) 
dt 



f(Xi(t),t) 



3 



CLijXj(t) - s{xi(t) - s(t)) 

i = 1,2,3 (23) 



Then, the coupling matrix becomes 

-10.0 5.0 0.1 
A = I 5.0 -11.0 6.0 
0.1 6.0 -6.1 

The largest eigenvalue Ai = —10.11 of the matrix A. There- 




Fig. 4. The solution of coupled system 1231 is pinned to s(t) with a single 
controller. 



IV. Conclusions 

Synchronization is an important research field in sciences 
and applications. How to pin a complex network to a specified 
solution (or an equilibrium point) of the uncoupled system is 
of great significance. However, in practice, the state variables 
of some nodes are not observable or measured. Therefore, 
we have to investigate the possibility of pinning a complex 
network by adding controllers to those nodes, which can be 
measured or controlled. In this paper, we prove rigorously 
that we can pin a complex network by adding a single linear 
controller to one node with symmetric or asymmetric coupling 
matrix. We also discuss how to pin the complex network with 
nonlinear coupling. Simulations also verify our theoretical 
results. 



Fig. 5. The solution of coupled system 1211 converges to s(t) through an 
asymmetric coupling 
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fore, Afe + cAi < for k = 1, ■ • • ,n, and the conditions of 
Theorem 2 are satisfied. The coupled system is synchronized 

, \ ~y 3 _ iM*)-s(t)n 

to sit). Figure 4 shows . The solution of 

coupled system (I23> is pinned to s(t) with a single controller 

B. The coupling matrix is asymmetric 

In this simulation, we take the following asymmetric matrix 




as the coupling matrix. As for the coupling strength, we pick 
e = 2. So 




Direct calculation indicates that the left eigenvector corre- 
sponding to eigenvalue is £ = (|, |, |) T and /ii = 0.0718, 
so we can choose c = 72, such that theorem 4 is satisfied and 
the simulation is shown by Figure 5. 



